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(57) Abstract: The invention relates to a 
method of determining impulse responses 
from a medium (2) in relation to the 
transmission of waves between different 
points (Tl-TN). The inventive method 
consists in: transmitting waves in the 
medium by generating, at each point i, 
signals ei(t) each comprising n elementary 
signals with respective frequencies which 
are spaced apart in pairs by an interval 5f 
and which are different from the frequencies 
of the elementary signals corresponding 
to the other points; receiving signals rj(t) 
at points j after the transmission of the 
aforementioned waves in the medium; and 
calculating each impulse response hij(t) 
from a correlation signal between signal 
ei(t) transmitted at point i and signal rj(t) 
received at point j. 

(57) Abrege : Proc6d6 pour determiner des 
reponses impulsionnelles d'un milieu (2) 
vis-a-vis de la transmission d'ondes entre 
diflferents points (Tl-TN), on emet des 
ondes dans le milieu en gSnerant en chaque 
point i des signaux ei(t) constitues chacun 
de n signaux elernentaires de frequences 
respectives espacees deux a deux d'un 
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pas 8f et distinctes des frequences des signaux 616mentaires correspondant aux autres points, puis on capte des signaux rj (t) en des 
points j apres transmission desdites ondes dans le milieu, et on calcule chaque reponse impulsionnelle hij(t) a partir d'un signal de 
correlation entre le signal ei(t) 6mis au point i et le signal rj(t) capte* au point j. 
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Procede pour determiner des reponses impulsi onnelles d'un 

milieu vis-a-vis de la transmission d' ondes entre 

dif f erents points . 

5 La presente invention est relative aux procedes 

pour determiner des reponses impulsionnelles d'un milieu 
vis-a-vis de la transmission d' ondes entre dif f erents 
points . 

Plus particulierement, 1' invention concerne un 
10 procede pour determiner des reponses impulsionnelles d'un 
milieu vis-a-vis de la transmission d' ondes entre 
dif f erents points, proced^ comprenant : 

(a) au moins une etape d' emission au cours de 
laquelle on emet des ondes dans le milieu en generant des 

15 signaux ei(t) a partir d'un nombre N de points d' emission 
appartenant au milieu, ou N est un nombre entier au moins 
egal a 2 et i est un iridice compris entre 1 et N qui 
designe l'un desdits N points d' emission, 

(b) au moins une etape de reception au cours de 
20 laquelle on capte des signaux rj (t) a partir desdites ondes 

apres transmission dans ledit milieu, en un nombre M de 
points de reception appartenant au milieu (ces points de 
reception peuvent etre confondus avec les points 
d' emission), ou M est un entier naturel non nul et j est un 
25 indice compris entre 1 et M qui designe l'un desdits M 
points de reception, 

(c) et au moins une etape de determination 
desdites reponses impulsionnelles hij (t) entre chaque point 
d' emission i et chaque point de reception j a partir des 

30 signaux <§mis ei(t) et captes rj (t) . 

Les reponses ' impulsionnelles peuvent etre 
determinees aux fins de focalisation d' ondes dans le 
milieu, par exemple aux fins d'imagerie du milieu ou de 
communication selective entre differents points du milieu, 

35 notamment dans le cadre de la mise en oeuvre de la 
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methode "D.O.R.T." (Decomposition de l'Operateur 

Retournement Temporel) decrite notamment par Prada et al. 
(C. Prada and M. Fink, "Eigenmodes of the time reversal 
operator: a solution to selective focusing in multiple- 
5 target media." Wave Motion, 20, pp 151-163 (1994)). 

Le processus d' acquisition des responses 
impulsionnelles entre les points en question du milieu est 
toutefois relativement lent, dans la mesure ou il est 
effectue sequentiellement , pour les differents points 

10 d' emission i les uns a la suite des autres, comme decrit 
par exemple dans le document WO-A-02/32 316 dans le cadre 
d'un procede de focalisation d' ondes acoustiques autre que 
la methode "D.O.R.T." susmentionnee. 

Cette lenteur allonge le temps de mise en ceuvre du 

15 procede, et lorsque le milieu est changeant, elle ne permet 
pas d'obtenir les reponses impulsionnelles assez vite pour 
pouvoir etre valablement utilisees ensuite dans ce milieu. 

La presente invention a notamment pour but de 
pallier ces inconvenients . 

20 A cet effet, selon 1' invention, un procede du genre 

en question est caracterise en ce qu'au cours de 1' etape 
(a), on fait §mettre simultanement les signaux ei(t) par 
lesdits N points d' emission, ces signaux ei(t) ayant une 
duree T et etant chacun une somme de n signaux elementaires 

25 sensiblement monochromatiques, de meme amplitude et de 
frequences respectives f 0 ,i+Jc.6f / ou f 0 ,i est une frequence 
predetermines propre au point i, k est un nombre entier 
compris entre 0 et n, n est un nombre entier au moins egal 
a 2 et 5f est un pas frequentiel predetermine, les 

30 frequences f 0 ,i respectives propres aux differents points i 
etant distinctes les unes des autres et comprises dans une 
bande de frequences de largeur 5f, 

et en ce qu'au cours de l'dtape (c) , on calcule chaque 
reponse impulsionnelle hij (t) a partir d'un signal de 
35 correlation entre le signal ei(t) emis au point i et le 
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signal rj (t) capte au point j. 

Grace a ces dispositions, on calcule simultanement 
les reponses impulsionnelles entre N points d' emission et M 
points de reception (eventuellement confondus avec les 
5 points d' emission), ce qui accelere le processus 
d' acquisition des reponses impulsionnelles. 

Dans divers modes de realisation du procede selon 
1' invention, on peut eventuellement avoir recours en outre 
a l'une et/ou a 1' autre des dispositions suivantes : 
10 - les frequences f 0 ,i respectives propres aux 

differents points i sont separees deux a deux par un ecart 
5f/N ; 

- au cours de l'etape (c) , ledit signal de 
correlation est fenetre au moyen d'une fonction porte n(t) 

15 de largeur l/5f ; 

- au cours de l'etape (c) , les reponses 
impulsionnelles hij (t) sont determinees par la formule : 

hij(t) = n(0- $eK# ~ f)Tj{&)dO ; 

- les ondes transmises dans le milieu entre les 
20 points d' emission et les points de reception sont des ondes 

acoustiques ; 

- au cours de l'etape (a), le milieu oil sont 
emises les ondes est reverberant ; 

- le pas frequentiel 5f est inferieur ou egal a 
25 1/t, ou x est la dispersion temporelle du milieu (c'est-a- 

dire la duree du signal capte apres emission d'un signal 
impulsionnel dans le milieu) ; 

- le pas frequentiel 5f est sensiblement egal a 
1/t, ou x est la dispersion temporelle du milieu ; 

30 - • la duree T est au moins egale a N/5f ; 

la duree T est au moins egale a N.x, ou t est la 
dispersion temporelle du milieu ; 

- les signaux elementaires presentent des phases 

aleatoires ; 
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- les ondes sont emises avec une certaine bande 
passante, les frequences fOi comprennent une frequence 
minimale fO et le nombre n est determine pour que la bande 
de frequence comprise entre fO et f 0+ [ (n+1) . 5f ] recouvre 
5 sensiblement ladite bande passante ; 

les points de reception sont confondus avec les 
points d' emission . 

D'autres caracteristiques et avantages de 
l 1 invention apparaitront au cours de la description 
10 suivante de deux de ses modes de realisation, donnes a 
titre d'exemples non limitatifs, en regard des dessins 
joints . 

Sur le dessin, la figure 1 est une vue schematique 
d'un dispositif adapte pour mettre en ceuvre la presente 
15 invention, dans un mode de realisation. 

Dans l f exemple non limitatif represents sur le 
dessin, le procede selon un mode de realisation de 
1' invention est mis en ceuvre par un dispositif 1 qui est 
adapts pour emettre et recevoir des ondes dans un milieu 2 
20 homogene ou non, generalement reverberant, qui peut etre 
solide, liquide ou gazeux. Les ondes en question sont de 
preference acoustiques et la description ci-apres ne fera 
reference qu' a de telles ondes acoustiques, mais lesdites 
ondes pourraient eventuellement etre d' une autre nature, 
25 par exemple de nature electromagnetique . 

Le dispositif 1 comporte par exemple un ordinateur 
3, ou tout autre dispositif de commande similaire, qui 
commande les emissions d' ondes acoustiques et leur ecoute 
dans le milieu 2. A cet effet, 1' ordinateur 3 peut par 
30 exemple commander un dispositif electronique dedie 4 qui 
peut notamment comporter une unite centrale electronique 
CPU associee a une memoire centrale MEM, cette unite 
centrale CPU commandant elle-meme N voies (N etant un 
entier au moins egal a 2) reliees chacune a un transducteur 
35 Tl-TN permettant d' emettre et de capter des signaux 
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acoustiques dans le milieu 2. Toutes les voies sont 
synchronisees entre elles. 

Chaque voie de mesure du dispositif 4 peut 
comporter par exemple un echantillonneur El -EN relie au 
5 transducteur correspondant Tl-TN, et chaque echantillonneur 
peut etre lui-meme relie a une memoire associee a une unite 
centrale, referencees conjointement Cl-CN. 

Le dispositif 1 peut etre utilise notamment aux 
fins d' imagerie du milieu 2, ou encore pour etablir une 
10 communication selective avec un point dudit milieu, ou 
autre, par exemple en utilisant la methode "D.O.R.T." 
(Decomposition de 1'Operateur Retournement Temporel) 
d6crite notamment par Prada et al. (C. Prada and M. Fink, 
"Eigenmodes of the time reversal operator: a solution to 
15 selective focusing in multiple-target media." Wave Motion, 
20, pp 151-163 (1994) ) . 

Dans tous les cas, il est necessaire de determiner 
les reponses impulsionnelles entre les transducteurs Tl-TN. 

A cet effet, l'ordinateur 3 commande le dispositif 
20 1 de fagon qu'il suive successivement les etapes 
suivantes : 

(a) au moins une etape d' emission au cours de 
laquelle on emet des ondes dans le milieu en g^nerant des 
signaux ei(t) a partir des N transducteurs Tl-TN (ces N 

25 transducteurs peuvent le cas echeant ne constituer qu'une 
partie des transducteurs relies au dispositif 4, auquel cas 
les reponses impulsionnelles relatives a 1" ensemble des 
transducteurs seraient acquises en plusieurs passes 
successives, ces passes restant toutefois en nombre tres 

30 inferieur au nombre total de transducteurs), i etant un 
indice compris entre 1 et N qui designe 1'un desdits N 
transducteurs ; 

(b) au moins une etape de reception au cours de 
laquelle les N transducteurs Tl-TN captent des signaux 

35 rj (t) a partir desdites ondes apres transmission dans ledit 
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milieu, ou j est un indice compris entre 1 et N designant 
l'un desdits transducteurs ; 

(c) et au moins une etape de determination 
desdites reponses impulsionnelles hij (t) entre chaque point 
5 d' emission i et chaque point de reception j a partir des 
signaux emis ei(t) et captes rj (t) . 

Pour permettre d'acquerir en une seule fois tous 
les signaux necessaires a la determination des reponses 
impulsionnelles concernant les N transducteurs 

10 susmentionnes (c ' est-a-dire tout ou partie des 
transducteurs reli§s au dispositif 4, comme explique ci- 
dessus), tous les signaux ei(t) sont emis simultan<§ment au 
cours de 1' etape (a) par les N transducteurs, mais ces 
signaux emis sont orthogonaux entre eux de fagon que les 

15 informations propres a chaque point d' emission puissent 
ensuite etre separees les unes des autres dans les signaux 
captes rj (t) . 

Ainsi, chaque signal emis ei(t) est egal a la somme 
de n signaux elementaires sensiblement mono ch r omati que s, de 
20 meme amplitude et de frequences respectives f 0 ,i+k.5f , ou : 

f 0/ i est une frequence predetermines propre au 

transducteur i, 

k est un nombre entier compris entre 0 et n, 
n est un nombre entier au moins egal a 2, 
25 - et 5f est un pas frequentiel predetermine, 

avantageusement inferieur ou egal a 1/x et de preference 
sensiblement egal a 1/t, ou t est une valeur moyenne de la 
dispersion temporelle du milieu 2, c' est-a-dire la duree du 
signal capte • apres emission d'un signal impulsionnel 
30 (fonction de Dirac) dans le milieu 2. 

Les frequences f 0 ,i respectives propres aux 
differents points i sont distinctes les unes des autres et 
comprises dans une bande de frequences de largeur 5f : 
avantageusement, ces frequences f 0/ i sont separees deux a 
35 deux par un ecart 6f/N. 
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De plus, si BP est la bande passante avec laquelle 
les ondes acoustiques sont emises et regues par les 
transducteurs, et si on nomme fO la plus faible des 
frequences fOi, alors fO et n peuvent avantageusement etre 
5 determines pour que la bande des frequences- des signaux 
elementaires (bande comprise entre fO et f 0+ [ (n+1 ) . 5f ] ) 
recouvre sensiblement ladite bande passante. 

Avantageusement, les signaux elementaires composant 
les signaux ei(t) presentent des phases aleatoires. 
10 Tous les signaux ei(t) ont une meme duree T, qui 

est avantageusement au moins egale a N/5f , et de preference 
au moins egale a N.t". 

Au cours de l'etape (c) , on calcule chaque reponse 
impulsionnelle hij (t) a partir d' un signal de correlation 
15 entre le signal ei(t) emis au point i et le signal rj (t) 
capte au point j . 

Dans l'exemple considere ici, ce signal de 

correlation vaut: 

cij(t)= $ei(0).rj(t-0)d0 (cette integrale est faite par 

20 exemple sur une duree au moins egale a T, dont le debut 

coincide avec le debut du signal rj (t) ) . 

Avantageusement, ce signal de correlation est 

fenetre au moyen d' une fonction porte n(t) de largeur l/6f 

(valant 1 pendant une periode de l/5f et 0 en dehors de 
25 cette periode), auquel cas les reponses impulsionnelles 

hij (t) peuvent avantageusement etre determinees par la 

formule : 

Ay(0 = n(0.ciy(0 . 

On notera que, dans l'exemple particulier 
30 represents sur le dessin et decrit ci-dessus, les points de 
reception j auxquels on capte l'onde acoustique sont 
confondus avec les points d' emission et constitues par les 
transducteurs Tl-TN, mais les points de reception en 
question peuvent le cas echeant etre constitues par un 
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second reseau de transducteurs (non represente), distinct 
du reseau de transducteurs Tl-TN. 

Dans le cas. general, les ondes acoustiques qui sont 
generees par I 1 emission des signaux ei(t) sont done captees 
5 en M points de reception, M etant un nombre entier non nul, 
lesquels M points de reception peuvent etre soit confondus 
avec les N points d 1 emission (auquel cas M = N) , soit non 
totalement confondus avec les points d 1 emission, soit 
entiereraent distincts des points d 1 emission. 
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RE VEND I CAT I ONS 

1. Procede pour determiner des reponses 
impulsionnelles d'un milieu (2) vis-a-vis de la 
5 transmission d' ondes entre differents points (Tl-TN) , 
procede comprenant : 

(a) au moins une etape d' emission au cours de 
laquelle on emet des ondes dans le milieu (2) en generant 
des signaux ei(t) a partir d'un nombre N de points 

10 d' emission (Tl-TN) appartenant au milieu, ou N est un 
nombre entier au moins egal a 2 et i est un indice compris 
entre 1 et N qui designe l'un desdits N po±nts d' emission, 

(b) au moins une etape de reception au cours de 
laquelle on capte des signaux rj (t) a partir desdites ondes 

15 apres transmission dans ledit milieu, en un nombre M de 
points de reception (Tl-TN) appartenant au milieu, ou M est 
un entier naturel non nul et j est un indice compris entre 
1 et M qui designe l'un desdits M points de reception, 

(c) et au moins une etape de determination 
20 desdites reponses impulsionnelles hij (t) entre chaque point 

d' emission i et chaque point de reception j a partir des 
signaux emis ei(t) et captes rj (t) , 

caracterise en ce qu'au cours de 1' etape (a), on fait 
emettre simultanement les signaux ei(t) par lesdits N 

25 points d' Emission (Tl-TN) , ces signaux ei(t) ayant une 
duree T et etant chacun une somme de n signaux elementaires 
sensiblement monochromatiques, de meme amplitude et de 
frequences respectives f 0 ,i+k.5f / ou f 0 ,± est une frequence 
predetermines propre au point i, k est un nombre entier 

30 compris entre 0 et n, n est un nombre entier au moins egal 
a 2 et 5f est un pas frequentiel predetermine, les 
frequences f 0 ,i respectives propres aux differents points i 
etant distinctes et comprises dans une bande de frequences 
de largeur 5f, 

35 et en ce qu' au cours de 1' etape (c) , on calcule chaque 
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reponse impulsionnelle hij (t) a partir d' un signal de 
correlation entre le signal 'ei(t) emis au point i et le 
signal rj (t) capte au point j . 

2. Procede selon la revendication 1, dans lequel 
5 les frequences fo,i respectives propres aux differents 

points i sont separees deux a deux par un ecart 5f/N. 

3. Procede selon la revendication 1 ou la 
revendication 2, dans lequel au cours de l'etape (c) , 
ledit signal de correlation est fenetre au moyen d'une 

10 fonction porte n(t) de largeur 1/Sf. 

4. Procede selon la revendication 3, dans lequel au 
cours de l'etape (c) , les reponses impulsionnelles hij (t) 
sont determinees par la formule : 

hijXO = n(f). \ei{9 - t).rj(G)de . 

15 5. Procede selon l'une quelconque des 

revendications precedentes, dans lequel les ondes 
transmises dans le milieu entre les points d' emission et 
les points de reception sont des ondes acoustiques. 

6. Procede selon l'une quelconque des 

20 revendications precedentes, dans lequel au cours de l'etape 
(a), le milieu ou sont emises les ondes est reverberant. 

1. Procede selon l'une quelconque des 

revendications precedentes, dans lequel le pas frequentiel 
5f est inferieur ou egal a ou t est la dispersion 

25 temporelle du milieu. 

8. Procede selon la revendication 7, dans lequel le 
pas frequentiel 5f est sensiblement egal a 1/t, ou t est la 
dispersion temporelle du milieu. 

9. Procede selon l'une quelconque des 
30 revendications precedentes, dans lequel la duree T est au 

moins egale a N/Sf . 

10. Procede selon l'une quelconque des 
revendications precedentes, dans lequel la duree T est au 
moins egale a N.t, oft x est la dispersion temporelle du 
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milieu. 

11. Procede selon 1 ' une quelconque des 

revendications precedentes, dans lequel les signaux 
elementaires presentent des phases aleatoires. 
5 12. Procede selon l'une quelconque des 

revendications precedentes, dans lequel les ondes sont 
emises avec une certaine bande passante, les frequences fOi 
comprennent une frequence minimale fO et le nonibre n est 
determine pour que la bande de frequence comprise entre fO 
10 et f0+ [ (n+1) -5f ] recouvre sensiblement ladite bande 
passante . 

13. Procede selon l'une quelconque des 
revendications precedentes, dans lequel les points de 
reception sont confondus avec les points d' emission (Tl- 
15 TN) . 
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